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5} ' ABSTRACT 

>• ; In this note we suggest that high- mass black holes; i.e., black holes of 

£^ ■ several solar masses, can be formed in binaries with low-mass main-sequence 

companions, provided that the hydrogen envelope of the massive star is removed 

^ in common envelope evolution which begins only after the massive star has 



finished He core burning. That is, the massive star is in the supergiant stage, 
which lasts only ~ 10 4 years, so effects of mass loss by He winds are small. 



Oh" Since the removal of the hydrogen envelope of the massive star occurs so late, 

it evolves essentially as a single star, rather than one in a binary. Thus, we can 
use evolutionary calculations of Woosley & Weaver (1995) of single stars. 

Using the Brown & Bethe (1994) upper limit of ~ 1.8 M Q for the 
(gravitational) compact core mass that can evolve into a low-mass black hole, 
^ ■ we find that high-mass black holes can be formed in the collapse of stars with 

ZAMS mass > 20 M . We somewhat arbitrarily take the upper limit for the 
evolution of the so-called transient sources to be ~ 35 M ZAMS mass. Mass 
loss by winds in stars sufficiently massive to undergo the LBV (luminous blue 
variable) stage may seriously affect the evolution of stars of ZAMS > 35 — 40 M Q , 
but we need calculations with improved mass loss rates before discussing these 
quantitatively. Both Portegies Zwart, Verbunt & Ergma (1997) and Ergma & 
van den Heuvel (1998) have suggested that roughly our chosen range of ZAMS 
masses must be responsible for the transient sources. We believe that the 
high-mass black hole limit of ZAMS mass ~ 40 M suggested by van den Heuvel 
& Habets (1984) and later revised to > 50 M (Kaper et al. 1995) applies to 
massive stars in binaries, which undergo RLOF (Roche Lobe Overflow) early in 
their evolution. We will not pursue this here because calculations with improved 
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He-star mass loss rates by wind, now being carried out, are necessary before 
quantitative results can be obtained. 

The most copious high-mass black holes of masses ~ 6 — 7 M Q have been 
found in the transient sources such as A0620. These have low-mass companions, 
predominantly of < 1 M , such as K- or M-stars. In the progenitor binaries the 
mass ratios must have been tiny, say q ~ 1/25. Normally such small ratios are 
thought to be rare; e.g. in binary evolution the companion distribution is often 
taken as dq, implying an very low probability of such a binary. 

In this note we follow the evolutionary scenario for the black hole binary 
of de Kool et al. (1987). We show that the reason for this small q- value lies 
in the common envelope evolution of the binary. The smaller the companion 
mass, the greater the radius R g the giant must reach before its envelope meets 
the companion. This results because the orbit of a low-mass companion must 
shrink by a large factor in order to expel the envelope of the giant, hence the 
orbit must initially have a large radius. (Its final radius must be just inside its 
Roche Lobe, which sets a limit to the gravitational energy it can furnish.) 

A large radius R g in turn means that the primary star must be in the 
supergiant stage. Thus it will have completed its He core burning while it is 
still "clothed" with hydrogen. This prevents excessive mass loss so that the 
primary retains essentially the full mass of its He core when it goes supernova. 
We believe this is why K and M-star companions of high-mass black holes are 
favored. 

We find that the black holes in transient sources can be formed from stars 
with ZAMS masses in the interval 20 — 35 M . The black hole mass is only 
slightly smaller than the He core mass, typically ~ 7 M Q . 

Subject headings: black hole physics — X-rays: bursts — stars: binaries 

1. DEVELOPMENT 

Our evolutionary scenario is essentially the same as that of de Kool et al. (1987) for 
the black hole binary A0620-00. We apply this scenario somewhat schematically for the 
range of ZAMS 20 — 35 M for the massive star. In Fig. [I] we show results of calculations 
of Fe core masses of single stars by Woosley & Weaver (1995). Somewhere around ZAMS 
mass 20 M Q the Fe core exceeds the Brown & Bethe (1994) limit of ~ 1.8 M Q for low-mass 
black-hole formation, so we take this as the beginning of high-mass black hole formation. 
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The detailed behavior of this curve should not be taken seriously, but the large increase 
around 20 — 25 M Q is of importance, as we discuss in detail later. 

We wish to carry out a population synthesis so that we can estimate the number of 
transient sources. We make roughly the same assumptions as Bethe & Brown (1998,1999), 
our massive star M B lying in mass somewhere in between the 0,B star progenitors of binary 
neutron stars and the progenitor of the massive black hole in Cyg X-l. Our low-mass 
companion is a main sequence star of mass ~ 1 M . Thus, the ratio 



is very small, and there will be great uncertainty in the initial number of binaries for such 
a small q ~ 1/25. We take the distribution as dq. The distribution in q is unknown for 
such low-mass companions as are involved here, but our results will show that the flat 
distribution in q is not unreasonable. We assume log a, where a is the semi-major axis of 
the orbit, to be uniformly distributed, over a total logarithmic interval of 7. Thus, the 
fraction of binaries in a given interval of In a is 



We take a the supernova rate to be ~ 0.02 per galaxy per year, somewhat larger than the 
0.015 given by Cappellaro et al. (1997) and somewhat smaller than the 0.025 assumed by 
Bethe & Brown (1998). The rate of supernovae (SN) is the same as the rate of birth of 
massive stars M > 10 M , so 



with n the Salpeter exponent which we take to be 1.5. We furthermore assume a binarity 
of 0.5. Once again, for our very small q values this is uncertain. 

We thus evolve as typical a ZAMS 25 M Q (star B) with companion ~ 1 M main 
sequence star (star A) as typical progenitor of the transient X-ray sources. The common 
envelope evolution can be done as in Bethe & Brown (1998). With M^. = 25 M Q and 
neglect of the accretion onto the main sequence mass M4, we find from Bethe & Brown 



d(f) 



d(\n a) 



(1.2) 
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(1.4) 



where Y = M B /a. Here the coefficient of dynamical friction Cd was taken to be 6, in 
the range of the 6-8 for supersonic flow with Mach number 3-10 (Ruffert 1994; Ruffert & 
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Arnett 1994). The result is relatively insensitive to q, the exponent 1.2 resulting from 

l + l/fo-l). 



Thus, in our case 
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We expect a ce ~ 0.5, under the assumption that the kinetic energy of the expelled envelope 
is equal to that it originally possessed in the massive star, but it could be smaller. From 
this we obtain 

* _ Mb^I - go(™. Me)™ 3 (16) 

where we have taken the He star mass Mb. to be 1/3 of Mg.. In order to survive spiral- in, 
the final separation af must be sufficient so that the main sequence star lies at or inside its 
Roche Lobe, about 0.2a/ if Ma = M . This sets a/ ~ 5i? = 3.5 x 10 11 cm and 

/0 ^\ ' 83 

ttl = 3.15 — xl0 13 cm. (1.7) 

\a ce J 

~ 2 A.U.. This exceeds the radius of the red giant tip in the more numerous lower mass 
stars in our interval, so the massive star must generally be in the supergiant phase where it 
meets the main sequence star, i.e., the massive star must be beyond He core burning. E.g., 
the red giant tip (before the He core burning) for a 20 M star is at 0.96 x 10 13 cm, for a 
25 M star, 2.5 x 10 13 cm (Schaller et al. 1992). These numbers are, however, somewhat 
uncertain. Notice that decreasing a ce will increase aj. Decreasing M A has little influence, 
because with the smaller stellar radius the minimum aj will decrease nearly proportionately. 
Note that neglect of accretion onto the main sequence star would change the exponent 0.83 
to unity, so accretion is unimportant except in increasing the final mass. 

Now a ZAMS 25 M star ends up at radius 6.7 x 10 13 cm (~ 2 a^) following He 
shell burning (Weaver, Zimmerman &Woosley 1978). Thus the interval between a« and 
6.7 x 10 13 cm is available for spiral-in without merger^] so that a fraction 
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6.7 



, ,, . 0.83 

3.15 



0.5 



0.11 (1.8) 



1 Note that envelope removal does not occur at the Roche Lobe on the thermal time scale r t h but at the 
low-mass star since the remaining lifetime of the giant is ~ 10 4 yrs, much shorter than r t h- 
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of the binaries survive spiral-in, but are close enough so that the main sequence star is 
encountered by the evolving H envelope of the massive star. The He core burning will be 
completed before the supergiant has moved out to ~ 2 A.U., so binaries which survive 
spiral-in will have He cores which burn as "clothed" , namely as in single stars. 



2. BIRTH RATE OF TRANSIENT SOURCES 

Given our assumptions in Section [I], the fraction of supernovas which arise from ZAMS 
stars between 20 and 35 M Q is 

0.20 (2.1) 



2 3/2 3.53/2 

where we have assumed the mass 10 M Q is necessary for a star to go supernova. A Salpeter 
function with index n = 1.5 is assumed here. Our assumption that the binary distribution 
is as dq is arbitrary, and gives us a factor 1/25 for a 1 M Q companion. Note, however, 
that had we included higher mass companions, the change in the final projected number 
of transient sources would not be order of magnitude, because the hydrogen burning time 
goes inversely as mass squared. Thus, for supernova rate 2 per century, our birth rate for 
transient sources in the Galaxy is 

2 x 10~ 2 x 0.5 x 0.11 x 0.20 x 0.04 ~ 8.8 x 10"V _1 (2.2) 



where 0.5 is the assumed binarity, 0.11 comes from eq. ( |1.8| ), and the final (most uncertain) 
factor 0.04 results from a distribution flat in q and an assumed 1 M companion star. 

In order to estimate the number of transient sources with black holes in the Galaxy, we 
should know the time that a main sequence star of mass ~ 1 M Q transfers mass to a more 
massive companion. For a main-sequence donor, the mass transfer rate is ~ 10~ 10 M Q yr _1 , 
almost independent of donor mass (Verbunt & van den Heuvel 1995). As mass is transferred, 
the mass of the donor decreases and with it the radius of the donor. Quite a few low-mass 
X-ray binaries have X-ray luminosities that imply accretion rates in excess of 10 -10 M yr _1 , 
leading to suggestions of additional mechanisms for loss of angular momentum from the 
binary, to increase mass transfer. Verbunt & Zwaan (1981) estimate that magnetic braking 
can boost the transfer of mass in a low-mass binary. We somewhat arbitrarily adopt an 
effective mass transfer rate of 10~ 9 yr -1 for main sequence stars and 10 _8 yr _1 for the two 
systems that have subgiant donors (V404 Cyg and XN Sco 94). In order to estimate the 
number of high-mass black hole, main sequence star binaries in the Galaxy we should 
multiply the birth rate (O) times the 10 9 yrs required, at the assumed mass loss rate, to 



strip the main sequence star, obtaining 8800 as our estimate. Not all of these will fill their 
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Roche Lobes. Those that do not may not now be visible, but will be later, as they begin 
evolving as one of those with subgiant donors. The fact that two of seven observed binaries 
are subgiants, although the lifetime of the latter is two orders of magnitude less than the 
main sequence lifetime, suggests that some fraction of our 8800 estimated binaries do not 
fill their Roche Lobes. From the observed black-hole transient sources Wijers (1996) arrives 
at 3000 low-mass black hole sources in the Galaxy, but regards this number as a lower limit. 
Beginning from this, which we regard as an observational estimate, we note that the two 
subgiants in Table 1 involve the more massive F,G and A stars. These indicate an ~ 100 
times greater population of unevolved main sequence stars in this range which lie quietly 
inside their Roche Lobes. Thus, including the quiescent binaries might give as many as 
(2/7) x 100 x 3000 or ~ 10 5 additional binaries, suggesting that our above estimate may 
be an order of magnitude too low. Estimates of the number of transient sources are very 
uncertain, but it is clear that there are orders of magnitude more of them than of the 
Cygnus X-l type objects with high-mass black hole and massive star companion. 

If we assume that ZAMS masses ~ 10 — 18 M evolve into a neutron star, we should 
have ~ 3 times more neutron stars than high- mass black holes (see eq. (|2.1| )). The upper 
limit follows from our belief that SN (1987A) with progenitor ~ 18 M Q ZAMS went into 
a low-mass black hole, following the scenario of Brown & Bethe (1994). On the basis of 
a Monte Carlo calculation using the kick velocities of Cordes & Chernoff (1997) we find 
that ~ 1/2 of the binaries containing He-star, low- mass main sequence companion (with 
M ~ 1 M ) will be disrupted in the explosion. Thus, we find only a slightly higher birth 
rate for LMXB's with neutron stars than with black holes, although the numbers could be 
equal to within our accuracy. The LMXB's with neutron stars tend to be much brighter 
than those with black holes, indicating an order of magnitude greater transfer rate. With 
the correspondingly shorter main-sequence lifetime, this would give us several hundred 
LMXBs with neutron stars, a factor of several greater than the observed number, ~ 130. 
Given the lifetime of ~ 10 8 years of a LMXB and accretion at roughly the Eddington rate 
of ~ 10 -9 — 10 -8 M yr _1 , it is reasonable that some neutron stars accrete a reasonable 
fraction of a solar mass M (van den Heuvel 1995). We expect the masses of these to 
exceed the Brown & Bethe (1994) limit of 1.5 M for maximum neutron star mass, and 
evolve into low-mass black holes. The low-mass black-hole, low-mass main sequence star 
systems would not be seen. 
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3. ESTIMATED MASSES OF THE BLACK HOLES IN TRANSIENT 

SOURCES 

As we showed below eq. (jl.7p , the He core of the massive star will in general be 
uncovered only after He core burning is completed. The remaining time for He burning (in 
a shell) will be short, e.g., for a 20 M© ZAMS star it is only 1.4 x 10 4 years (Schaller et al. 
1992). Therefore the mass loss by wind after uncovering the He core will not be large, and 
when the star finally becomes a supernova, its mass will be almost equal to the He core of 
the original star. The latter can be calculated from 

M He ~ 0.08(M ZAMS ) L4 (3.1) 

so for ZAMS masses 20 — 35 M Q M# e will lie in the interval ~ 5.3 — 11.6 M Q . In fact, the 
lower limit looks a bit small, because the He core of the ~ 18 M Q progenitor of 1987A is 
generally taken as ~ 6 M . 

Bailyn et al. (1998) find the black hole masses in transient sources to be clustered 
about ~ 7 M Q , except for V404 Cyg which has a higher mass. This is in general agreement 
with our scenario, because most of the black holes will come from the more numerous stars 
of ZAMS mass not far from our lower limit of ~ 20 M . Two points are important to note: 

1. Not much mass can have been lost by wind. Naked He stars have rapid wind loss. 
However in our scenario the He star is made naked only during He shell burning and 
therefore does not have much time (< 10 4 yrs) to lose mass by wind. 

2. There are good reasons to believe that the initial He core will be rotating (cf. 
Mineshige, Nomoto & Shigeyama 1993). The way in which the initial angular 
momentum affects the accretion process has been studied by Mineshige et al. (1997) 
for black hole accretion in supernovae. In general accretion discs which are optically 
thick and advection dominated are formed. The disc is hot and the produced energy 
and photons are advected inward rather than being radiated away. The disc material 
accretes into the black hole at a rate of > 10 e MEdd for the first several tens of days. 
Angular momentum is advected outwards. Our results show that little mass is lost, 
because the final ~ 7 M Q black hole masses are not much less massive than the 
He core masses of the progenitors, and some mass is lost by wind before the core 
collapses. The latter loss will not, however, be great, because there is not much time 
from the removal of the He envelope until the collapse. 

Accretion of the He into the black hole will differ quantitatively from the above, but we 
believe it will be qualitatively similar. The fact that the helium must be advected inwards 
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and that little mass is lost as the angular momentum is advected outwards is extremely 
important to establish. This is because angular momentum, essentially centrifugal force, 
has been suggested by Chevalier (1996) to hold up hypercritical accretion onto neutron 
stars in common envelope evolution. (Chevalier (1993) had first proposed the hypercritical 
accretion during this evolutionary phase to turn the neutron stars into black holes, the 
work followed up by Brown (1995) and Bethe & Brown (1998).) However, once matter is 
advected onto a neutron star, temperatures > 1 MeV are reached so that neutrinos can 
carry off the energy. The accreted matter simply adds to the neutron star mass, evolving 
into an equilibrium configuration. Thus, this accretion does not differ essentially from that 
into a black hole. In either case of neutron star or black hole an accretion disc or accretion 
shock, depending on amount of angular momentum, but both of radius ~ 10 11 cm, is first 
formed, giving essentially the same boundary condition for the hypercritical accretion in 
either case, black hole or neutron star. Thus, the masses of the black holes in transient 
sources argue against substantial inhibition of hypercritical accretion by jets, one of the 
Chevalier (1996) suggestions. 

Measured mass functions, which give a lower limit on the black hole mass are given in 
Table |] Only GRO J0422+32 and 4U 1543-47 have a measured mass function < 3 M . 
Results of Callanan et al. (1996) indicate that the angle % between the orbital plane and the 
plane of the sky for GRO J0422+32 is i < 45°, and recent analysis by Orosz et al. (1998) 
indicate that the angle i for 4U 1543-47 is 20° < i < 40°. So both GRO J0422+32 and 4U 
1543-47 also contain high- mass black holes. 

4. GENERAL DISCUSSION 

There is agreement (Portegies Zwart et al. 1997; Ergma & van den Heuvel, 1998) that 
in order to make enough transient sources the progenitors of the black holes must begin at 
relatively low masses, ~ ZAMS 20 M Q . We take the upper limit somewhat arbitrarily to 
be ~ 35 M about where rapid mass loss occurs and stars may enter into the LBV phase. 
Our upper limit is relatively unimportant since most of the stars considered will lie near the 
lower limit. 

Based on the observations of Kaper et al. (1995) that the companion is a hypergiant, 
Ergma & van den Heuvel (1998) argue that the progenitor of the neutron star in 4U1223-62 
must have a ZAMS mass > 50 M . Brown, Weingartner & Wijers (1996), by similar 
argumentation, arrived at ~ 45 M , but then had the difficulty that 4U1 700-37, which 
they suggested to contain a low-mass black hole appeared to evolve from a lower mass star 
than the neutron star in 1223. Wellstein & Langer (1999) suggest the alternative that in 
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1223 the mass occurs in the main sequence phase (Case A mass transfer) which would be 
expected to be quasi conservative. They find that the progenitor of the neutron star in 1223 
could come from a mass as low as 26 M . This is in agreement with Brown et al. (1996) 
for conservative mass transfer (their Table 1), but these authors discarded this possibility, 
considering only RLOF (Case B mass transfer) in which case considerable mass would be 
lost. 

Wellstein & Langer (1999) are in agreement with Brown et al. (1996) that 4U1700-37 
should come from a quite massive progenitor. Conservative evolution here is not possible 
because of the short period of 3.4 days (Rubin et al. 1996). The compact object mass is 
here 1.8 ± 0.4 M (Heap & Corcoran 1992). Brown et al. (1996) suggest that the compact 
object is a low-mass black hole. The upper mass limit for these was found by Brown & 
Bethe (1994) to be ~ 1.8 M , as compared with an upper limit for neutron star masses of 
~ 1.5 M . Thus, there seems to be evidence for some ZAMS masses of ~ 40 — 50 M ending 
up as low-mass compact objects, whereas we found that lower mass stars in the interval 
from ~ 20 — 35 M ended up as high-mass black holes. In this sense we agree with Ergma 
& van den Heuvel (1998) that low-mass compact object formation "is connected with other 
stellar parameters than the initial stellar mass alone." We suggest, however, following 
Brown et al. (1996) that stars in binaries evolve differently from single stars because of 
the different evolution of the He core in binaries resulting from RLOF in their evolution. 
Namely, "naked" He cores evolve to smaller final compact objects than "clothed" ones. 

In fact, this different evolution of binaries was found by Timmes, Woosley & Weaver 
(1996). They pointed out that stars denuded of their hydrogen envelope in early RLOF 
in binaries would explode as Type lb supernovae. They found the resulting remnant 
gravitational mass following explosion to be in the interval of 1.2 — 1.4 M , whereas in 
exploding stars of all masses with hydrogen envelope (Type II supernova explosion) they 
found a peak at about 1.28 M , chiefly from stars of low masses and another peak at 
1.73 M more from massive stars. Our Fe core masses in Fig. [I| come from essentially the 
same calculations, but the "Remnant" masses of Woosley & Weaver (1995) are somewhat 
greater than those used by Timmes et al. (1996). In fact, the differences between the masses 
we plot and those of Timmes et al. come in the region ~ 1.7 — 1.8 M (gravitational). This 
is just in the Brown & Bethe (1994) range for low- mass black holes. It may be that some 
of the stars with low-mass companions evolve into low-mass black holes. Presumably these 
would give lower luminosities than the high-mass black holes, although at upper end of the 
mass range we discuss 4U1700-37 seems to be an example of such a system. Of course here 
the high luminosity results from the high mass loss rate of the giant companion. There are 
substantial ambiguities in fallback, etc., from the explosion. Our point in this paper is that 
most of the higher mass single stars 20 — 35 M go into high mass black holes. (The Brown 
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& Bethe (1994) limit for low-mass black hole formation is ~ 1.5 — 1.8 M gravitational, but 
there is some give and take in both lower and upper limit. Also the stars are not all the 
same. In particular different metallicities will give different wind losses.) 

In determining the upper mass limit for which a low-mass compact object can result 
from binary evolution with RLOF (Case B mass transfer) the Brown et al. (1996) scenario 
must be revised because it is now realized that He-star wind loss rates employed by Woosley, 
Langer & Weaver (1995) were too high. With lower rates the differences in behavior 
between "clothed" " and "naked" He core evolutions will be diminished because much of the 
difference arises from the mass losses and how it affects the convective 12 C burning as we 
discuss below. 

In the mass range ~ 20 — 35 M , the compact objects resulting from naked He stars 
(Type lb SN explosions) are sufficiently far below the maximum neutron star mass that they 
will still remain neutron stars when better (lower) He star mass loss rates are employed. 
Indeed, for our 25 M star, Wellstein and Langer (1999) find that halving the mass loss 
rate lowers the mass of the CO core by only 2.5 %. However, at about the upper end of this 
range of masses, the outcome may, indeed, be changed. We will somewhat arbitrarily focus 
on the mass region ~ 40 M , studying how lower mass loss rates might affect the outcome. 

Until recently W.-R. wind loss rates were taken from observed winds which originated 
chiefly from free-free scattering. These depend quadratically on density. Because 
of "dumpiness" in the winds, the mass loss rate was overestimated. Polarization 
measurements of the Thomson scattering, which depend linearly on the wind, give 
substantially lower mass rates, in approximate agreement with the rates that would be 
deduced from the observed rate of increase in orbital periods for spherical mass loss 

M IP , . 

H = T- 

In V444 Cygni P = 0.202 ± 0.018 s yr" 1 was obtained by Khaliullin et al. (1984) and 
M WR = 9.3 ± 0.5 M by Marchenko et al. (1994), resulting in 

M dyn = 1.03 x 10~ 5 M^r- 1 . (4.2) 

This is to be compared with the 

M = 0.75 x 10" 5 Mqjy- 1 (4.3) 

obtained by St. -Louis et al. (1993). In later work Moffat & Robert (1994) arrive at a mean 
of (0.7 ± 0.1) x 10~ 5 M yr _1 . The mass loss rate employed by Woosley, Langer & Weaver 
(1995) was that of Langer (1989). Specifically, M = -kM 2 - 5 , with M in M and M in 
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M yr _1 and k = 6 x 10~ 8 as long as the carbon surface mass fraction does not exceed 0.02 
and k = 10~ 7 afterwards. Choosing an average of k = 8 x 10~ 8 we find the WLW rate to 
be 2.1 x 10~ 5 M yr _1 , a factor of 2 larger than Mdyn- Given the many uncertainties in our 
estimate, we feel the range of two to three times less than the Langer (1989) mass loss rate 
to be reasonable. It should also be remembered that stars vary substantially in metallicity, 
and that there will be a range of variation even of those in the Galactic disc. Of course 
stars in the metal poor Magellanic clouds should have substantially lower winds. 

In the region of ZAMS ~ 40 M calculation of Wellstein & Langer (1999) show that 
both the final He core mass and C/O core mass increase ~ 23% in the mass loss rate is 
halved. Thus, reducing the mass loss rate by a factor of 2 does not increase the final mass 
by a similar factor. The main reason is that if the mass loss is reduced, then the He stars 
remain somewhat more massive, thus also more luminous and therefore have a higher mass 
loss rate than had their mass been reduced earlier. There are also other, less important, 
feedbacks. 

In the case of ~ 60 M stars, the He core and C/O cores increase ~ 31% if the mass 
loss rate is halved and a factor of ~ 1.8 — 1.9 if it is cut by 1/4. In fact, in this case, a 
reasonable interpolation formula for these lower mass loss rates is 

TT = (l-33) 1/2/ (4.4) 

where M was calculated with the Langer (1989) rates and / is the fractional decrease in 
winds from these rates. Wellstein & Langer (1999) also evolve a 36 M He core, which 
would come from an ~ 85 M ZAMS star, with / = 1/4 of the Langer et al (1989) wind 
loss rates. Their final 7.5 M He star mass is not enough for Cyg X-l. However, Woosley, 
Langer & Weaver (1993) with the Langer 1989 mass loss rate found a 9.71 M He core 
in evolving a ZAMS 85 M star with mass loss. In the latter case the WNL phase lasted 
> 1/3 of the W.-R. phase; i.e., the He core initially had some hydrogen envelope. This 
illustrates the order of uncertainties that may come in the evolution of very massive stars, 
arising from mass loss in the LBV phase. 

Although the final He core and C/O core scale by roughly the same factor the Fe 
core, crucial for the compact object masses, is not expected to do the same. Firstly, rather 
trivially, even if the baryon number Fe core did scale in the same way, one would expect 
corrections downwards in core mass from binding energy correction. These decreases in 
gravitational mass would be greater for higher mass cores, so the Fe cores would scale 
somewhat less rapidly than the C/O cores. 

More important for the mass of the Fe core is the ZAMS mass at which the convective 
carbon burning is skipped, because, as seen in the single stars in Fig. [1], a big jump in Fe 



-12- 



core masses occurs here. 

The convective carbon burning phase (when it occurs) is extremely important in 
presupernova evolution, because this is the first phase in which a large amount of entropy 
can be carried off in uu-pa.ii emission, especially because this phase is of long duration. 
The reaction in which carbon burns is 12 C(a,7) 16 (other reactions like 12 C+ 12 C would 
require excessive temperatures). The cross section of 12 C(a, 7) 16 is still not accurately 
determined; the lower this cross section the higher the temperature of the 12 C burning, and 
therefore the more intense the vv emission. With the relatively low 12 C(a,7) 16 ratesQ 
determined both directly from nuclear reactions and from nucleosynthesis by Weaver & 
Woosley (1993), the entropy carried off during 12 C burning in the stars of ZAMS mass 
< 18 M is substantial. The result is rather low-mass Fe cores for these stars, which can 
evolve into neutron stars. Note that in the literature earlier than Weaver & Woosley (1993) 
often large 12 C(a,7) 16 rates were used, so that the 12 C was converted into oxygen and 
the convective burning did not have time to be effective. Thus its role was not widely 
appreciated. 

Of particular importance is the ZAMS mass at which the convective carbon burning 
is skipped. In the Woosley & Weaver (1995) calculations this occurs for single stars at 
ZAMS mass 19 M Q but with a slightly lower 12 C(a,~f) l6 rate it might come at 20 M Q or 
higher (Brown 1997). As the progenitor mass increases, it follows from general polytropic 
arguments that the entropy at a given burning stage increases. At the higher entropies 
of the more massive stars the density at which burning occurs is lower, because the 
temperature is almost fixed for a given fuel. Lower densities decrease the rate of the triple-a 
process which produces 12 C relative to the two-body 12 C (a, 7) 16 which produces oxygen. 
Therefore, at the higher entropies in the more massive stars the ratio of 12 C to 16 at the 
end of He burning is lower. The star skips the long convective carbon burning and goes 
on to the much shorter oxygen burning. Oxygen burning goes via 16 0+ 16 giving various 
products, at very much higher temperature than 12 C(a,7) 16 and much faster. Since 
neutrino cooling during the long carbon-burning phase gets rid of a lot of entropy of the 
core, skipping the convective carbon burning phase leaves the core entropy higher and the 
final Chandrasekhar core fatter. 

In Fig. [1] the large jump in compact object mass in single stars at ZAMS mass ~ 19 M Q 
is clearly seen. From our discussion in the last section we see that this is just about at the 



2 Weaver & Woosley use S(E) = 5(300) = 170 keV barns, remarkably close to the 169 ± 55 keV barns 
arrived at by Barnes (1995). Given the large uncertainty (stemming chiefly from that in the E2 rate, the 
good agreement may be somewhat accidental. 
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point where our Fe core mass goes above ~ 1.8 M Q and, therefore, above this mass one 
would expect single stars to go into high-mass black holes. Arguments have been given that 
SN (1987A) with progenitor ZAMS mass of ~ 18 M evolved into a low-mass black hole 
(Brown & Bethe 1994). We believe from our above arguments and Fig. |I| that soon above 
the ZAMS mass of ~ 18 M Q , single stars go into high-mass black holes without return 
of matter to the Galaxy. Thus, the region of masses for low-mass black hole formation 
in single stars is narrow. The precise upper mass limit is not clear, but certainly in the 
range of > 20 M ZAMS. This is in agreement with Wellstein & Langer (1999) who find a 
minimum black hole progenitor mass for single stars of 21 M . 

Thus far our discussion has been chiefly about single stars, in which the He burns while 
clothed by a hydrogen envelope. 

In binary evolution, if the hydrogen envelope is removed before the helium core burning 
in either Case A or Case B (RLOF) mass transfer, the resulting "naked" He star burns 
quite differently from a "clothed" one. Weaver & Woosley (1993) find that the convective 
carbon burning tends to be skipped when the central 12 C abundance at the end of helium 
core burning is less than ~ 15%. In Woosley, Langer & Weaver (1995) with the Langer 
(1989) helium mass loss rate, this central 12 C abundance is 34% for a ZAMS 40 M Q star. 
With mass loss rate decreased to half, it is still 33% (Wellstein & Langer, 1999) hardly 
changed. However, as noted earlier, the CO mass for a ZAMS 40 M star is increased from 
2.33 M to 2.87 M Q with halved mass loss rate. Provisionally Wellstein & Langer (1999) 
have suggested that the magnitude of the CO mass chiefly determines the fate of the star, 
and that the 2.87 M CO core could go into a black hole. However, halving the mass loss 
rate of the naked He star decreases the central carbon abundance at the end of He core 
burning hardly at all. Thus, there will still be a long period of 12 C convective core burning. 
As seen from Fig. [I] for less massive stars in the mass region ~ 20 M , this can easily lower 
the compact core mass by ~ 0.5 M . On the other hand, for each 1 M added to the 
He envelope at the envelope at the time of SN explosion, an additional ~ 3 x 10 49 ergs is 
needed to expel it in either a prompt or delayed supernova explosion. This is ~ 2% of the 
SN explosion energy, so would work towards formation of a high-mass black hole. 

Given the increase in He and CO cores from the decreased mass loss rates, we estimate 
that stars of ZAMS masses ~ 40 M which lose their H envelopes by wind will end up as 
low-mass black holes, 1700-37 being an example. The possible mass ranges in which this 
can happen cannot be determined until the dynamic evolution of the CO cores formed with 
the lowered He-star wind rates of Wellstein & Langer (1999) is carried out. 

It seems clear that with the lower metallicity in the LMC and consequently lower mass 
loss rates that the 12 C convective core burning will not be skipped so the Fe cores in the 
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> 40 M region of ZAMS masses will be larger, and will evolve into black holes. This may 
help to explain why there are two HMXB's with high-mass black holes, LMC X-l and LMC 
X-3, in the LMC, whereas Cyg X-l is the only clear example in the disc. 

5. CONCLUSION 

We have shown that it is likely that single stars in the range of ZAMS masses 
~ 20 — 35 M© evolve into high-mass black holes without return of matter to the Galaxy. 
This results because at mass ~ 20 M the convective carbon burning is skipped and this 
leads to substantially more massive Fe cores. Even with more realistic reduced mass loss 
rates on He stars, it is unlikely that stars in this mass range in binaries evolve into high-mass 
black holes, because the progenitor of the compact object when stripped of its hydrogen 
envelope in either Case A (during main sequence) or Case B (RLOF) mass transfer will 
burn as a "naked" He star, ending up as an Fe core which is not sufficiently massive to form 
a high-mass black hole. 

In the region of ZAMS mass ~ 40 M , depending sensitively on the rate of He-star 
wind loss, the fate of the primary in a binary may be a low-mass black hole. We are unable 
to pin down the limit for high-mass black hole formation until better mass loss rates are 
determined. 

In our estimates we have assumed the Brown & Bethe (1994) estimates of 1.5 M for 
maximum neutron star mass and 1.5 — 1.8 M for the range in which low- mass black holes 
can result. 

In our evolution of the transient sources using Case C (during He shell burning) mass 
transfer, almost the entire He core will collapse into a high-mass black hole, explaining 
the more or less common black hole mass of ~ 7 M for these objects, with the possible 
exception of V404 Cygni where the mass may be greater. Our evolution gives an explanation 
for the seemingly large gap in black- hole masses, between the > 1.5 M for the black hole 
we believe was formed in 1987A and the ~ 1.8 M black hole we suggest in 1700-37 and the 
~ 7 M in the transient sources. 

We note that following the removal of the H envelope by Case C mass transfer, the 
collapse inwards of the He envelope into the developing black hole offers the Collapsar 
scenario for the most energetic gamma ray bursters of Woosley (1993) and MacFadyen & 
Woosley (1999). Especially if our higher estimate of ~ 10 5 of the high-mass black-hole, 
main-sequence star binaries is roughly correct, a simple estimate (Brown, Lee & Wijers 
1999) shows that they would be the largest population of possible progenitors. 
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Table 1: Parameters of suspected black hole binaries with measured mass functions (Wijers 
1996, Chen et al. 1997, Bailyn et al. 1998, Orosz et al. 1998, Bailyn, private communication). 
N means nova, XN means X-ray nova. Numbers in parenthesis indicate errors in the last 
digits. 
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(*) A much higher inclination for A0620 has been claimed by Haswell et al. (1993) of up 
to i=70. In this case, the lower limits on the component masses would be M x > 3.8 and 
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Fig. 1. — Comparison of the compact core (Remnant) masses resulting from the evolution 
of single stars (filled symbols), case of solar metallicity of Woosley & Weaver (1995). 



